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16448 TWO BED-LOAD FORMULAS: EVALUATION 


KEY WORDS: Bed load; Beds under water; Channel beds; Channels 
(waterways); Comparative studies; Formulas (mathematics); Instrumentation; 
River beds; Sands; Sediments; Sediment transport 


ABSTRACT: evaluation formulas presented for the prediction bed load 
sand-bottomed channels. The evaluation done measuring both the hydraulic data 
needed for the formulas and the actual sediment transport rate, and comparing the 
predicted with the observed values. The results indicate that the Toffaleti procedure 
predicts the transport rate with sufficient accuracy useful for most engineering 
purposes. 


REFERENCE: Amin, Magdy (Lect., Dept. Civ. Engrg., Ain Shams Univ., Cairo, 
Egypt), and Murphy, Peter J., “Two Bed-Load Formulas: Evaluation,” Journal 
the Hydraulics Division, ASCE, Vol. 107, No. HY8, Proc. Paper 16448, August, 1981, 
pp. 961-972 


16422 PROFILES STRATIFIED FLOWS 


KEY WORDS: Density; Fluid mechanics; Open channel flow; Saline water; 
Stratification; Turbulent flow; Velocity; Velocity measurement 


ABSTRACT: The Monin-Obukhov analysis applied stably stratified, two- 
dimensional, hydraulically rough, turbulent, fully developed, open-channel flows. The 
buoyancy flux expressed function the velocity and density gradients for cases 
which the buoyancy flux may vary over the depth. Resulting vertical velocity profile 
expressions were verified with experimental data collected wide laboratory flume 
with salt-stratified flows. Sidewall effects stratified flows are significant narrow 
flumes and are magnified highly stratified flows, noted review previous 
works. 


REFERENCE: McCutcheon, Steve (Hydrologist, U.S. Geological Survey, Gulf 
Coast Hydroscience Center, NSTL Station, Miss., 39529), Velocity Profiles 
Stratified Flows,” Journal the Hydraulics Division, ASCE, Vol. 107, No. HY8, 
Proc. Paper 16422, August, 1981, pp. 973-988 


16449 MODEL FOR FLOW OVER DREDGED CHANNELS 


KEY WORDS: Approach channels; Channels (waterways); Computation; 
Computer models; Dredging; Eddies; Estuaries; Flow patterns; Mathematical 
models; Nonuniform flow; Turbulence 


ABSTRACT: The mathematical model nearly horizontal flows applied for the 
study the nonuniform flow over steep-sided, dredged, port approach channel. 
mixed finite-difference/finite-element method used for the solution. The eddy 
viscosity parameter computed means mixing length and turbulence 
kinetic energy turbulence closures. The comparison the computed flow magnitudes 
(velocity, eddy viscosity, shear stresses, and turbulence kinetic energy) the 
laboratory measured ones demonstrate the predictive character the model. 


REFERENCE: Koutitas, (Lect., Hydraulics Dept., School Tech., Thessaloniki 
Univ., Thessaloniki, Greece), and O’Connor, B., “Turbulence Model for Flow over 
Dredged Channels,” Journal the Hydraulics Division, ASCE, Vol. 107, No. HY8, 
Proc. Paper 16449, August, 1981, pp. 989-1002 
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16454 ROUTING DIVERGING-CONVERGING WATERSHEDS 


KEY WORDS: Conversion; Diversion; Floods; Mathematical models; 
Overland flow; Rainfall; Rainfall-runoff relationships; Surface dynamics; 
Surface runoff; Watersheds; Wave equations 


ABSTRACT: For diverging-converging surface, the kinematic wave equations are 
solved using implicit finite difference method. 
representation applicable natural watershed that has leaf-shaped geometry. 
The results are compared with other idealistic representations. These surface, 
diverging surface, and converging surface. The diverging-converging model gives 
hydrograph values between those converging and diverging models, depending 
the geometric parameters. general, the geometry watershed has important 
effects the runoff hydrograph for both equilibrium and partial equilibrium 
conditions. 


REFERENCE: Agiralioglu, Necati (Asst. Prof. Civ. Engrg., Technical Univ. 
Istanbul, Istanbul, Turkey), “Water Routing Diverging-Converging Watersheds,” 
Journal the Hydraulics Division, ASCE, Vol. 107, No. HY8, Proc. Paper 16454, 
August, 1981, pp. 1003-1017 
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Magdy Amin,' and Peter Murphy,’ ASCE 


INTRODUCTION 


There are many formulas for the prediction the rate sediment transport 
rivers. However, the predictions not agree with one another, and 
they not define the transport rate satisfactory way. This problem was 
considered Vanoni (9) and Shen (7), and state-of-the-art recommendations 
were made. The focus this paper the evaluation two these predictive 
formulas when they are applied sand-bottomed channels. The formulas selected 
for study were chosen both for their perceived importance and correspond 
the limitations the measurement procedure. The first writer (1) has reported 
the evaluation equations. The best two the 10, the Meyer-Peter and 
Miiller formula (5) and the Toffaleti formula (8), will treated herein. 

Most the sediment discharge formulas predict bed-load discharge, because 
bed load logically dependent the local properties the channel. The 
prediction suspended sediment transport more difficult since the transport 
any cross section depends the balance between the particles going into 
suspension and the deposition from suspension over long upstream reach. 
Therefore, bed-load formulas were chosen. 

The evaluation done measuring both the hydraulic data needed for the 
formulas and also the actual sediment transport rate, and then comparing the 
predicted with the observed values. The bed-load formulas use the local mean 
velocity, flow rate, bed shear stress, stream power predict the transport 
rate particular site river. Since the equations assume uniform steady 
flow, apply them river channel necessary find: (1) reach 
the river that roughly straight and has fairly uniform width and depth; 
and (2) period when the flow steady. Even under these conditions the 
wide variation the depth with the distance across cross section makes 
desirable divide the cross section into subsections, usually done 
when current-meter measurements were made. The test site was chosen meet 
these conditions, and the cross section was divided into five subsections. 

Dept. Civ. Engrg., Ain Shams Univ., Cairo, Egypt. 

Prof., Dept. Civ. Engrg., Univ. Massachusetts, Amherst, Mass. 01003. 

Note.—Discussion open until January 1982. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication July 22, 1980. 
This paper part the Journal the Hydraulics Division, Proceedings the American 


Society Civil Engineers, Vol. 107, No. HY8, August, 1981. ISSN 0044- 
796X /0008-0961 $01.00. 
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Test 


The sediment discharge data were taken Fish Creek, near Rome, N.Y. 
using the compartmented sediment trap described detail Ref. The test 
site shown Fig. was located one side island that was 
being formed the center wide crossover the meandering channel. 
The mean water-surface slope, S,,, was measured with point gages installed 
the channel, 343 apart, cross sections nearly equal area. The creek 
bed was covered with combination ripples and dunes the typical manner 
sand-bottomed channels. The wavelength and amplitude the dunes were 
roughly and 0.3 respectively, while the wavelength and amplitude 
the ripples were approx 0.3 and 0.03 The average value the sieve 
diameter, the sand the site was 0.4 mm. 


VELOCITY 
PROFILE 


VERTICALS 


CROSS SEDIMENT (6) 
SECTIONS 


ISLAND 


SECTIONS 
FIG. Creek Field Site 


Detailed flow measurements were made two cross sections the channel, 
shown Fig. Velocity profiles were taken the verticals shown, 
and samples the moving sediment were collected the base each vertical. 
The velocity was measured with Price current meter and the sediment was 
caught compartmented sediment traps. The grain size the bottom and 
the samples was measured with standard sieves. The water temperature was 
measured with mercury thermometer, and the duration the sediment collection 
period was observed with stopwatch. 

The placement the sediment traps the stream bed presented some difficulty. 
first, digging tower similar offshore oil drilling platform was used. 
Its operation was slow that only one trap could used each day. Next, 
digging frames were used shallow water where wading was possible, and 
this procedure was sufficiently rapid that trap samples with their corre- 
sponding velocity profiles could obtained day. Finally, the digging frames 
were used divers the cold, deep, fast-flowing water the spring runoff. 
That process was roughly half fast the wading procedure. Due the 


lo 
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presence suspended sediment, few samples were taken these higher flow 
rates. 
The operation the traps and the collection the auxiliary data required 


four people. outline the procedure presented herein. More detail 
available Ref. 


Install the digging frame the channel bed near the middle dune. 

Fill the trap with water and place the closed trap the frame. 

Wait until the natural state the ripples returns the sand bed above 
the trap. 

Measure the velocity profile, U(z), and depth, each trap location. 

Measure the mean channel slope, S,,, and water temperature, 

Open the trap’s cover and collect sediment sample during period 
sufficient duration gather 0.5 kg. During this period observe the trap periodically 
insure proper operation. The collection period varied from 6.5 

Close the trap and remove from the stream bed without spilling the 
sample. the trap was open for more than day, repeat the auxiliary data. 


Dry, weigh, sieve, and reweigh the samples using standard grain size analysis 
techniques. 


The field data were analyzed order calculate the hydraulic variables 
typically used describe river. The friction velocity, V,, and the height 
zero velocity, z,, the logarithmic velocity profile were obtained from 
linear regression analysis the velocity profile data for each trap location. 
The correlation coefficients were usually above 0.90 for these analyses. The 
mean velocity, was obtained from the numerical integration velocity profile 
over each trap’s vertical. The local bed shear stress, was calculated from 
and the Darcy-Weisback friction-factor, was calculated from 
The velocity profile data were checked for accuracy calculating the von 
Karmen constant, from the mean velocity equation 


The mean value was 0.404 with standard deviation 0.035, indicating 
that the measured velocity distributions were the expected logarithmic profiles. 
The sediment samples collected the containers each trap were weighed 
and the spacial distribution the sediment the trap was checked for accuracy 
calculating the jump length probability density function outlined Ref. 
Samples with inconsistent jump length patterns were discarded. The mass 
transport rate, g,, was calculated from which total mass collected, 
container width; and sample collection time. After sieving the collected 
sand obtain its grain size distribution, the mass components, with container 
number, and grain size, were also calculated and checked for consistency. 


Test 


The data taken over 2-yr period are presented Table where the transport 
rate calculated from each individual sample presented, along with the corre- 
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TABLE 1.—Experimental Data 


centi- centi- 
centi- per per 
meters second 


milli- 
meters 


964 HY8 
qs: in 
grams 
(1) (5) (6) (7) (9) 
100 3.4 4.5 16.4 1.7 0.33 
119 4.9 16.2 17.6 1.3 6.31 
1.9 3.6 21.4 0.5 0.28 0.012 
1.7 1.8 21.8 1.0 0.36 0.0054 
2.3 4.0 21.8 1.0 0.28 0.0030 
2.2 2.3 21.8 1.0 0.28 0.0062 
1.9 0.3 21.8 1.0 0.36 0.031 
2.6 4.1 23.8 0.9 0.30 0.010 
3.1 13.9 23.8 0.9 0.28 0.0052 
1.6 23.8 0.9 0.31 0.0029 
2.6 2.2 19.4 0.30 0.12 


BED-LOAD FORMULAS 
TABLE 1.—Continued 


sponding velocity profile parameters and other auxiliary data. Only the measured 
sediment transport rate given the table, although the fractional contributions 
the grain sizes the trap compartments were recorded. 

Three different reasons led the conclusion that the sediment transport 
data represented primarily bed load. First, visual observation the sand moving 
along the bottom showed sign suspension since the water was clear. 
Second, the spacial distribution the sand the trap indicated that the jump 
length only small fraction (~0.2%) the sample was greater than 1,000 
particle diameters. Finally, the ratio the friction velocity the fall velocity, 
was compared Bagnold’s suspension threshold (2), and 
was generally found less than that threshold. Nine samples failed meet 
the latter two criteria and were omitted from the data set. The resulting transport 
rate data are therefore considered bed-load. 

The quality the data set was studied further examining the power law 
relationship between the measured transport rate and the three hydraulic parame- 
ters, flow rate per unit width, stream power, and mean velocity. The resulting 
power laws (in units) were 


and 


The high exponents these power laws greatly amplify the effects errors 
the hydraulic parameters and Since the sediment samples were 
not much larger than the mass ripple and the von Karmen constant calculation 
indicated velocity error level 10%, smoothing was used provide larger 
samples and more accurate velocities for the evaluation the predictive 


2.0 0.4 23.2 0.34 0.011 
2.4 2.0 23.2 0.30 0.0039 
2.1 3.6 17.6 0.29 0.0020 
2.8 5.0 17.6 0.30 0.011 
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equations. The smoothing was done ordering the trap samples accordance 
with their sediment transport rates, and averaging each set five neighboring 
values these ordered sediment transport rates and the corresponding hydraulic 
parameters. The best the three smoothed relationships, that the mean 
velocity, shown Fig. The averaging improved the correlation coefficient 
this power law from 0.92-0.98. 

The comparison between the present data and previous bed-load measurements 
best done dimensionless format. The bed-load function, and the shear 


MEAN 
wu 


SEDIMENT TRANSPORT RATE-q, 


FIG. Transport Rate 


PLANE BED 


” 
” 
” 
z 
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DIMENSIONLESS TRANSPORT RATE— 


FIG. with Previous Data from Vanoni (9) 


stress function, are the principal dimensionless variables used describing 
bed load 


and [(y, — ¥) dso] 


Fig. the smoothed data are compared previous bed-load data, described 


- 
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Ref. for the movement coarser material over plane beds. The difference 
between the two data sets basically due Reynolds number effect. The 
lower Reynolds number the moving sand was shown 
Liu (3) permit the formation dunes and ripples, whereas the beds 
the earlier tests were planar. The shear stress applied the sand particles 
reduced the pressure drag these wave forms. Thus, this Reynolds 
number effect requires the separation the bed’s drag force into shear drag 


and form drag and causes the increased scatter seen the sand transport data 
Fig. 


AND 


The comparison between the sediment movement predicted the bed-load 
equations and the observed movement will presented ratio, that is, 
the predicted transport rate will divided the observed transport rate. 
Since additional details about these equations are available the works 
the original authors and Ref. only brief outline the prediction procedures 
will given herein. 

The Meyer-Peter and Miiller formula (5) for wide rectangular channel 


3/2 1/3 


which the specific weight water; the part the discharge 


apportioned the bed and considered responsible for the bed load; the 
total discharge; the total roughness coefficient for the bed; the 
flat-bed grain roughness; the water depth; the slope the energy 
the effective diameter the bed material; the acceleration due gravity; 
and the bed load weighed under water. 

The Meyer-Peter and Miiller equation the threshold type. The bed load 
proportional the difference between the mean shear stress acting the 
grain and the threshold shear stress, raised the 1.5 power. This threshold 
stress corresponds roughly Shield’s critical tractive stress. the application 
this equation the region near the threshold, the transport may appear 
negative. such cases, the authors the equation suggested replacing 
the coefficient 0.047 with 0.030 Eq. and that procedure was followed 
herein. 

The Meyer-Peter and Miiller equation was based data from experiments 
with steady, uniform flow flumes. The range sediment properties these 
tests was such that there was little suspended load. The effective diameter 
the sediment particles varied between 0.4 and mm. This range led 
flows both with and without bed waves and the bedform correction ratio, 
shown the formula varied between 0.5 and 1.0. This wide range 
sediment diameters also covered the cases both fully-rough bed flow 
and flow with viscous sublayer. the latter case, the authors the equation 
recommended the use the channel Reynolds number and the relative roughness 
the bed material find the friction factor, from the Nikuradse diagram. 
When the roughness, based this was combined the ratio, 
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this procedure was considered valid for the wide range sediment movements 
studied. 

The extensive variation the data base the Meyer-Peter and Miiller equation 
permits the application that equation the sediment movement Fish Creek. 
The equation was used predict the sediment movement based the smoothed 
hydraulic data for the creek’s flow. These predictions are compared with the 
measured sediment transport Fig. The geometric mean value the ratio 
predicted transport that observed was 1.92, showing the average tendency 
the formula predict excess transport. The ratio’s logarithmic standard 
deviation 0.29 shows that the scatter the agreement factor 2.0. 

The Toffaleti method for the calculation bed load (8) modification 
the Einstein procedure for the calculation both bed load and suspended 
load (3). focuses channels with sand bottoms, and assumes that part 
the sand moving bed load and part suspension. Thus the bed load, 
BL, found evaluating the mass flux, UC, 2D, and assuming that 


ERVED 


SEDIMENT TRANSPORT RATE— 


FIG. Meyer-Peter and Prediction with Observed Transport 
Rate 


this flux constant the region, 2D, near the bed 


which power law approximation the velocity profile; 
and the sediment concentration distribution. This velocity profile 
roughly 1/7 power law description based the turbulent boundary layer 
profile and adjusted for water temperature. The concentration particles 
grain size, determined equating the suspended load the region, 
with sediment load, GF, based Einstein’s bed-load function. 
The variable calculated from 


0.60| 


DTA 


the grain shear stress. The ratio, based Toffaleti’s unification 
Einstein’s and into the parameter, order more simply relate 


3 ° ° 6% 
q ike) 
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the grain shear stress parameter, the overall shear stress variable, 

The three principal differences between Toffaleti’s procedure and the Einstein 
method were: (1) The change the equation for the velocity distribution from 
the logarithmic velocity profile the power law profile; (2) the use instead 
and and (3) the use the suspended load the region 
R/11.24 evaluating the bed load. The changes were made order 
simplify the original procedure and optimize the agreement with the more 
extensive data base used Toffaleti (8). That base was composed 239 
points from six rivers and 282 points from five flume studies. Like the data 
base the Meyer-Peter and Miiller formula, this data base provided broad 
range hydraulic variables and sediment transport rates. 

The Toffaleti method was applied the smoothed Fish Creek data, and 
the resulting ratio the predicted observed bed load shown Fig. 
The geometric mean value that ratio was 1.15, showing close agreement 
between the Toffaleti predictions and the measured values. The ratio’s logarithmic 
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SEDIMENT TRANSPORT RATE— 
FIG. Toffaleti Prediction with Observed Transport Rate 


standard deviation 0.39 shows scatter this agreement factor 


The goal this paper, the evaluation bed-load equations, requires accurate 
data for both the channel’s hydraulic parameters and the sediment transport 
rate. The experimental procedure presented here met that need measuring 
the flow and sediment transport detail, using variety internal accuracy 
checks, and averaging the data reduce its scatter. The best evidence 
the quality the resulting data the correlation coefficient 0.98 between 
the bed load and the mean velocity. 

The comparisons between predicted and measured bed loads show that the 
Toffaleti method the better predictor the bed load sand-bottomed channels. 
However, obvious from Fig. that the scatter the agreement between 
the predicted and observed values much higher than the scatter 
the basic data, shown Fig. Since both the Toffaleti method and the 
best power law description the current data are based principally the 
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mean velocity, the additional scatter could associated with the effects 
the grain size, water depth, and temperature, and the slope the local energy 
grade line used the Toffaleti procedure, with the randomness the bed 
load. The Meyer-Peter and Miiller predictions were significantly higher than 
the measured sediment transport rates. Further, their method predicted zero 
movement many cases where movement was observed. This effect the 
threshold type predictive equation especially troublesome when bed forms 
are present, since the separation form drag from bed particle drag makes 
the location the threshold more difficult. These problems demonstrate that 
the Meyer-Peter and Miiller equation should not used for sand bottomed 
channels. 
The conclusions this evaluation are: 


The measurement procedure, using the velocity profiles and the compart- 
mented sediment trap, provides the accuracy and the self-consistency needed 
for reliable bed-load measurement, provided that sufficient samples are taken 
permit averaging. 

The Toffaleti bed-load procedure predicted the mean bed load the 
sand-bottomed channel Fish Creek with sufficient accuracy useful for 
most engineering purposes. 


Bed-load equations the threshold type, like the Meyer-Peter and 
equation, are not useful for sand-bed channels. 
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The following symbols are used this paper: 


BED-LOAD FORMULAS 


bed roughness parameter (8); 
width compartment trap; 

bed load (8); 

sediment concentration (8); 

velocity coefficient (8); 

sediment concentration coefficient (8); 
grain size (8); 

effective diameter bed material (5); 
grain sieve diameter; 

constant integration (8); 

friction factor; 

sediment load based Einstein’s (8); 
bed load weighed under water (5); 
acceleration gravity; 

water depth (5); 

depth open channel flow; 

number container; 

number grain size; 

flat-bed grain roughness (5); 

mass trap; 

mass grain size, compartment, 
total water discharge (5); 

water discharge responsible for bed load (5); 
sediment’s mass transport rate per unit width; 
hydraulic radius (8); 

particle Reynolds number; 

slope energy grade line (5); 

mean slope channel; 

water temperature; 

duration interval when trap open; 
local velocity (8); 

mean velocity (8); 

velocity fluid height, above bed; 
mean velocity fiuid; 

friction velocity; 

fall velocity; 

fall velocity particle; 

elevation above bed (8); 

elevation point (8); 

exponent velocity distribution (8); 
exponent concentration distribution (8); 
roughness length velocity profile; 
specific weight sediment and water; 
von Karmen’s constant; 

kinematic viscosity; 
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density fluid; 

dimensionless shear stress; 

bed shear stress; and 

dimensionless sediment transport rate. 


972 HY8 

= 

= 


AUGUST 1981 


VERTICAL VELOCITY PROFILES 
STRATIFIED 


INTRODUCTION 


Flows stratified heat, dissolved solids, and suspended solids are quite 
common the natural environment. These include flows the atmospheric 
boundary layer, rivers, lakes, estuaries, and oceans. The vertical transport 
the stratifying agent and momentum depend the shape the vertical velocity 
profile. 

The study described this paper focuses the shape the vertical velocity 
profile unidirectional, 2-dimensional, turbulent, fully developed, continuously 
stratified flow. For this case, the turbulence the flow generated the 
movement fluid over rough boundary. The dimensional analysis Monin 
and Obukhov (10) employed extend the Prandtl-von Karman velocity law 
stratified flows. The buoyancy flux related the velocity and density 
gradient such that the original Monin-Obukhov similarity analysis extended 
flows which the buoyancy flux may vary over the depth. Emphasis 
placed open-channel flows; therefore, flows laboratory flume stratified 
salt are used for verification. 


Previous Work 


Several studies have considered turbulent stratified flows over rough 
boundary, but the results these studies are conflicting. define and resolve 
this conflict, useful introduce schematic definition expected vertical 
velocity profiles stratified flows, shown Fig. Stable stratification describes 
flow which lighter layers fluid move over heavier layers and vice versa 
for unstable stratification. Neutral stratification implies well mixed homoge- 
neous flow. 

Fig. follows from extensive studies the 2-dimensional, stratified, atmo- 
spheric boundary layer (18). Schiller and Sayre (14), Hwang (6), the writer 
(8), Odd and Rodger (12), Vanoni (20), and Einstein and Chien (3) have measured 

U.S. Geological Survey, Gulf Coast Hydroscience Center, NSTL Station, 
Miss. 


Note.—Discussion open until January 1982. extend the closing data one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication April 1980. 
This paper part the Journal the Hydraulics Division, Proceedings the American 
Society Civil Engineers, Vol. 107, No. HY8, August, 1981. ISSN 0044- 
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velocity profiles stably stratified wide open-channel flows that resemble curve 
type Fig. Odd and Rodger measured vertical velocity profiles wide 
estuary during ebb tide. Schiller and Sayre and Hwang studied heat stratified 
flows flume; the writer studied salt stratified flow flume; and Vanoni 
and Einstein and Chien studied suspended-sediment stratified flows flumes. 

Flume studies French (4) and Majewski (9,10) narrow flumes not 
agree with the studies previously mentioned. French studied stable salt-stratified 
flows and Majewski studied stable heat-stratified flows. French and Majewski 
measured velocity profiles stable stratified flows that have the general shape 
curve rather than curve 

This conflict can resolved considering the width-to-depth ratios and 
the level stratification. From the summary the various studies stratified 
flows given Table evident that nonrepresentative velocity profile shapes 
occur open channels with width-to-depth ratios 4.3 and less. Studies 
Majewski are enlightening. Majewski ran control experiments with unstratified 
flows that document sidewall influence and strong secondary currents his 
narrow channels. hypothesized that the nonrepresentative velocity profiles 


3 


| UNSTRATIFIED 
2 STABLY STRATIFIED 
3 UNSTABLY STRATIFIED 


VELOCITY VELOCITY 


FIG. 1.—Three Types Velocity Profiles Expected 2-Dimensional Stratified Flows 


narrow stratified flows are result the three-dimensional effects similar 
those recorded Majewski’s control experiments. 

The work Schiller and Sayre also indicates that the level stratification 
may also lead nonrepresentative velocity profiles wider channels, shown 
Fig. Fig. shows the shape vertical velocity and temperature profiles 
along the channel center line from the work Schiller and Sayre. When 
high temperature difference exists over the depth, the velocity profile shapes 
follow those French and Majewski. Discounting entrance effects, 
hypothesized that stratification magnifies sidewall influence. 

The results Kantha al. (7) better illustrate the influence geometry 
stratified flows. Kantha al. studied stratified flow annular tank. 
Turbulence was input through the water surface grid moving over the 
surface. The growth well-mixed layer below the surface was related 
the shearing rate. Fig. shows the results Kantha al. There are geometry 
effects gross Richardson number 36. higher Richardson numbers, 
sidewall effects the annular tank become significant. Therefore, geometry 
and the level stratification are important producing 2-dimensional, turbulent, 
stratified flows. 
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number expressions have been put forth describe vertical profiles 
stratified flows. The Monin-Obukhov similarity analysis (11) results velocity 
profile expression that has been verified Webb (21) and others (11) for the 
atmospheric boundary layer which the buoyancy flux essentially constant 
with distance above the bottom boundary. Hwang (6) applied the Monin-Obukhov 
analysis open-channel flows with very limited results. Hwang assumed 
constant buoyancy flux and theorized that the momentum transfer unaffected 
stratification. Schiller and Sayre (14) derived profile expression from the 
Reynolds equations that based the Boussinesq approximation (small density 
changes have inertia effects), and this expression predicts decrease 


TABLE 1.—Summary Studies Stably Stratified Flows 


Width-to- 
depth ratio 
(2) 


Stratifying 
agent 
(3) 


Experimental 
device 


(4) 


Investigator 


(1) 


French (4) 1.8-0.8 salt flume 
Majewski (9) heat hydraulic model 
Majewski (10) 1.24, 2.07,3.10 heat flume 
Hino, al. (5) salt and sugar flume 


Schiller and Sayre 
(14) 

Hwang (6) 

Monin and Yaglom 
(11) 

Webb (21) 


5.0, 10.0 
5.0 


oo 


heat 
heat 
heat 


flume 
flume 
atmospheric bound- 
ary layer 
atmospheric bound- 
ary layer 


heat 


Chuang and Cermak 
(2) 


SethuRaman and 


heat wind 


Cermak (15) 9.0-16.0° heat wind tunnel 
Plate and Lin (13) heat wind tunnel 
Odd and Roger (12) 20-100 salt estuary 
Vanoni (19) 11.8-5.0 suspended flume 


sediment 
salt 


McCutcheon (8) 4.4, 8.0 


the average width the trapezoid sections. 
information. 


“Based the boundary layer thickness. 


flume 


velocities near the boundary due stable stratification. This prediction not 
borne out experimental evidence. 

French (4) also applied the Monin-Obukhov analysis develop velocity 
profile for stratified flows, but unlike Hwang (6) avoided the assumptions 
constant buoyancy flux and unaffected momentum transfer stratification. 
The development the following paragraphs follows French’s development 
many details, but French empirically modified the definition the Monin- 
Obukhov scaling length for stratified flows match his nonrepresentative narrow 


Curve 
(5) 
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VELOCITY, IN FEET PER SECOND 


~—— LOGARITHMIC VELOCITY PROFILE 
—— STRATIFIED FLOW VELOCITY PROFILE 
DIMENSIONLESS DISTANCE DOWNSTREAM 


DENSIMETRIC FROUDE NUMBER 
DEPTH=0.076m 


FIG. 2.—Velocity and Temperature Profiles Measured Schiller and Sayre (14) 
Flume Center Line 


TANK 


NORMALIZED ENTRAINMENT RATE 


BOTTOM WALL THE 


RATIO 


FIG. 3.—Entrainment Rate Function Depth-to-Width Ratio for Various Gross 
Richardson Numbers, for Stratification Entrainment Rate 
Normalized Values Depth/Width 0.5 and 1.41 cm/s. Adapted from 
Kantha al. (7) 
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flume velocity profiles. Therefore, the following analytical development strives 
rectify problems encountered French and Hwang. 


The vertical velocity profile unstratified flows has been the focus several 
studies (19,13), from which can concluded that the velocity profile 
hydraulically rough, 2-dimensional, unstratified flows can described the 
Prandtl-von Karman velocity law 


which bottom shear velocity; von Karman’s turbulence constant 
0.4); distance above bottom boundary with the bottom 
boundary; bottom roughness scale; and u(z) temporal average horizontal 
velocity some distance, from the bottom boundary. Using the Prandtl-von 
Karman velocity law basis, the similarity analysis Monin and Obukhov 
(11) employed herein develop vertical velocity profile expressions that 
are valid for open channel flows. 

Monin and Yaglom (11) note that turbulent, two-dimensional, stratified flow 
can characterized five parameters that include u,, and 


which kinematic viscosity and buoyancy flux, parameter dynamic 
significance defined 


which local acceleration gravity; depth-averaged density; and 
space- time-averaged fluctuations density and vertical velocity. 

Not all these five parameters are important determining the shape 
the vertical velocity profile. The kinematic viscosity not important turbulent, 
hydraulically rough flows, except very near the bottom boundary. The bottom 
roughness scale determines the magnitude velocity point but does not 
affect the shape the profile. Therefore, the shape the vertical velocity 
profile may described B). Based these parameters the Buckingham 
theorem predicts that single, dimensionless group can formed 


3 3- 


von Karman’s constant was introduced Obukhov (11, 475) and retained 
for consistency. 

This length scale, (the Monin-Obukhov scaling length), related the 
height dynamic sublayer next the bottom boundary which stratification 
has effect the dynamic structure the flow. this lowest layer the 
turbulence the flow overrides the potential energy effects stratification, 
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and the velocity profile can described Eq. Since this dimensional 
analysis has focused the shape the velocity profile, can 
hypothesized that dimensionless combination du/dz, u,, and related 


which required that when stratification disappears. This forces 


For small values can expanded power series 


which are numerical coefficients. Lacking sufficient information evaluate 
all the coefficients (21), truncated after the linear term, i.e., 


which constant for small values The accepted value from 
atmospheric studies 5.2 (20,11,18). Therefore, the velocity gradient 
can expressed 


pu, 


proceed, necessary evaluate the buoyancy term, These turbulent 
fluctuations can related the velocity and density gradients with the Prandtl 
mixing length hypotheses after giving some consideration the relationship 
between mass and momentum mixing lengths, 

general, function z/L. Lacking knowledge this 
relationship, becomes necessary consider small values which 
1.0 This compatible with earlier assumptions 
for that are valid for small The buoyancy term can then expressed 


buoyancy flux terms and the velocity and density gradient allows the 
scaling length, expressed depth-dependent variable. The Monin- 
Obukhov analysis can then extended stratified flows which the buoyancy 
flux varies over the depth. 

Eqs. and are then combined and algebraically manipulated solve for 
This results the following expression for velocity gradient 
stratified flow 
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which ca,. Thus the shape the velocity profile related the 
shape the density profile. Furthermore, for unstratified flow, which 
Eq. reduces the Prandtl-von Karman velocity law, Eq. 
can solved for velocity profile expression when the shape 
the density profile defined. Since many open-channel density profiles are 


nonlinear, useful consider general case which f(z), resulting 


which f(z) negative for stable density gradient, forcing the denominator 
less than one. Therefore, the right side Eq. greater than 


indicating that stable density gradient leads increase velocity over 
what would expected density stratification did not exist. Therefore, Eq. 
follows the expectations for stratified flows first expressed Taylor (16). 
Taylor found that stable stratification inhibits momentum transfer. 

the density profile can approximated straight line, 
the velocity profile can expressed 


(:) 


k? 2 


pu, 
important note that the denominator Eq. differs sign with 
the similar velocity profile expression French (4) because French incorrectly 
defined the Monin-Obukhov scaling length. The coefficient also different 
the factor French’s velocity profile expression indicates that stable 
stratification increases momentum transport. This must attributed sidewall 
influence. 

Having improved upon the approach Hwang and corrected sign difference 
introduced French, necessary verify Eqs. and with two-dimen- 
sional, stratified, open-channel data. The next section describes data collected 
the writer for that purpose. 


verify Eqs. and 14, series laboratory experiments was carried 
out flume 0.6 wide, 0.3 deep, and 18.3 long. Saline water was 
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injected underneath fresh water simultaneously with the use splitter plate 
produce stably stratified flows, shown Fig. The flume was operated 


TAILGATE 
POINT GAGES 


PITOT TUBE- 
CONDUCTIVITY 
OBE 


GATE VALVE 


FIG. 4.—Experimental Equipment Utilized for Stratified Flow Experiments 


FLOW DIRECTION 


BOTTOM ROUGHNESS CROSS SECTION 


0.035m 


FLOW 


UME BOTTOM 
SECTION 


FIG. Roughness Used Studies 


once-through mode with experimental times min and min. These 
times were determined discharge and saline water storage capacity. Experi- 
ments varied depth flow, discharge, bottom roughness patterns, and density 
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difference over the depth. Depths were approximately 0.08 and 0.14 with 
corresponding discharges 0.144 and 0.0074 Bottom roughness 
was two types. One consisted expanded metal mesh plaster lath 
and the other consisted 0.013-m wooden strips centers, shown 
Fig. Density profile shapes were controlled varying the initial density 
difference and controlling the interfacial mixing the end the splitter 
plate. Salinity was the sole stratifying agent; heat stratification was avoided. 
Density profiles varied from two-layer profiles linear density profiles. Values 
z/L varied from 0-4.5. 

Point velocities were measured over the depth the flume center line with 
Prandtl type 1/8 in. (3.2 mm) diameter pitot-static tube connected 
airless water column P90D Celesco variable reluctance, in. (25.4 mm) 
range, differential pressure transducer. (The use the brand name this report 
for identification purposes only and does not imply endorsement the U.S. 
Geological Survey.) Densities were deduced from saline solution electrical 
conductivity measured the same time with custom made probe described 
Anderson al. (1). The probe consisted hypodermic syringe one 
pole with insulated platinum wire the center extending out the needle 
make second pole. The change voltage across these two poles was 
related solution conductivity through extensive calibration. These profiles 
were simultaneously measured 7.0 downstream the splitter plate. 

Water surface slopes were computed from water-surface elevations measured 
with standard point gages [0.001 (0.000305 graduations] over 3.35-m 
reach approximately centered about the point where velocity and density profiles 
are measured. The gages were mounted over the flume center line. Discharge 
was measured with orifice meters the fresh-water and saline-water supply 
lines. 

Several observations give assurance that experiments were properly designed 
verify Eqs. and 14. Each experimental combination depth, discharge, 
and bottom roughness was applied first unstratified flow and then 
stratified flows. These unstratified control experiments, along with boundary 
layer calculations, indicated that the bottom boundary layer was fully developed 
for unstratified flows. Then velocity profiles, measured across the flume 
the test section, showed negligible influence from secondary circulation 
unstratified control experiments. Finally, fluorescent dye added the saline- 
water supply quickly became completely mixed over the depth unstratified 
control experiments. For stratified flows, the dye penetrated the water surface 
but remained stratified. Eddies generated the bottom were visually observed 
penetrate the less dense upper layer until eddies reached the water surface. 
The fluorescent dye and the density difference between the eddy and surrounding 
waters made possible track large eddies. 

The density difference between the upper layer and lower layer was decreased 
until eddies could observed moving from the bottom the water surface 
point well upstream the test section. This practice gave some assurance 
that velocity profiles were controlled bottom-generated turbulence. The 
concept uniform flow cannot applied stratified flows which the 
density profile decays the downstream direction. However, areas where 
bottom generated turbulence migrates over the depth flow causing density changes 
all depths, Eqs. and should valid. 
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Discharge, Temper- 
Depth, ature, Roughness 
Experiment cubic meters degrees scale, 
number per second meters Celsius meters 


(1) (2) (3) (6) 


(a) Wooden Strip 
RUN 2/13/3 
RUN 2/22/1 
RUN 2/13/2 
RUN 2/16/1 
RUN 2/15/1 
RUN 2/15/2 
RUN 2/20/3 
RUN 2/20/2 
RUN 2/21/2 
RUN 2/21/1 
RUN 2/20/4 


(6) Expanded Metal 
RUN 2/7/2 
RUN 2/6/1 
RUN 2/6/2 
RUN 2/5/3 
RUN 2/5/2 
RUN 2/5/1 
RUN 12/19/3 
RUN 1/22/2 
RUN 12/20/1 
RUN 11/10/2 
RUN 12/14/2 


Resutts 


Data collected for each experiment included velocity profiles and density 
profiles the flume center line 7.0 downstream the splitter plate, water 
surface slope over 3.4 reach approximately centered about that point, 
discharge, water temperature, and depth flow. Table presents these data 
and the resulting reduced data. 

The reduced data include the roughness scale, shear velocity, average velocity 
and density the vertical the flume center line, Reynolds number, and 
friction factor. The roughness scale and shear velocity were calculated using 
the two four lowest points the measured velocity profile for stratified 
flows. The roughness scale and shear velocity were also determined from 
unstratified control experiments for the purpose comparison. Shear velocities 
varied slightly with stratification, but control experiments served rough 
benchmarks for the rejection experiments with questionable shear velocity 
measurements. 


HY8 VERTICAL VELOCITY PROFILES 
draulic Data 


Shear Average Average 
velocity, velocity, density, Reynolds 
meters meters kilograms per number, 
per second per second cubic meter UD/v 
(7) (8) (9) (10) (11) 
Roughness 


999.83 
1,004.91 
1,004.40 

999.83 
1,003.29 
1,003.74 
1,001.35 
1,002.78 
1,002.78 
1,000.57 
1,001.42 


The roughness scale and shear velocity determinations should not affected 
stratification. This was checked comparing the velocity predictions 
Eq. (unstratified) and Eq. (stratified) all the points near the bottom 
used calculate the roughness scale and shear velocity. The greatest difference 
between and was found 4%. This indicates that stratification 
did not significantly bias the determinations the roughness scale and shear 
velocity. 

Twenty-two velocity profiles (109 velocity measurements) measured un- 
stratified and stratified flows are presented Ref. (8). Fig. represents 
typical example. Here measured velocity profile (circles) plotted semi- 
logarithmic and arithmetic scales and compared the predictions Eqs. 
14, and 12. This example indicates that Eqs. and are better predictors 
velocity profiles stratified flows than Eq. 

This better illustrated considering pairs dimensionless velocities 
consisting measured value and predicted value. Table gives the mean 
differences between measured dimensionless velocities and predicted dimension- 


0.0549 0.288 999.65 6.97 0.291 
0.0667 0.318 1,005.47 7.92 0.351 
0.0589 0.329 1,004.02 7.79 0.257 
0.0430 0.196 999.90 4.13 0.380 
0.0404 0.204 1,002.54 4.25 0.314 
0.0400 0.184 1,004.03 3.73 0.382 
0.0310 0.210 1,000.04 8.07 0.169 
0.0275 0.202 1,002.58 7.86 0.148 
0.0170 0.108 999.86 4.12 0.193 
0.0153 0.098 1,001.66 3.86 0.222 
0.0204 0.113 1,001.06 4.16 0.196 
Mesh Roughness 
0.0630 0.303 0.342 
0.0600 0.308 6.85 0.344 
0.0619 0.307 6.67 0.325 
0.0390 0.172 3.69 0.409 
0.0416 0.173 3.70 0.405 
0.0389 0.173 3.73 0.462 
0.0283 0.168 7.29 0.227 
0.0271 0.195 7.85 0.154 
0.0284 0.165 7.43 0.236 
0.0160 0.097 5.00 0.218 
0.0155 0.112 4.94 0.154 
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less velocities using: (1) Eq. (2) Eq. 14; and (3) Eq. 12. the 
means indicated significant difference measured velocity and predicted 
velocity, using Eq. whereas there not significant difference measured 
velocity and predicted velocity using Eq. Eq. 14. 

Prior evaluating Eqs. and 14, the Monin-Obukhov power series coefficient, 
was evaluated 5.8 0.3 from the laboratory data generated this study. 
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FIG. 6.—Typical Density Profile with Linear Approximation and Comparison Typical 
Measured Velocity Profile with Eqs. and Semilogarithmic and Arithmetic 
Scales 


The value the maximum value z/L, and the fit Eqs. and 
the measured data are interrelated and inseparable. Eqs. and can fit 
data with large (>0.0025) decreasing but this violates two assumptions 
Eqs. and 14: (i) Eq. limited small and (2) the variable 
Eq. not For these laboratory data, large can accommodated 
only varies with the maximum Therefore, the maximum was 
decreased incrementally 0.0025 until did not vary function 
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was evaluated changing that coefficient until the parameter 


was nearest zero. was the number velocity measurements. For 


max 


0.025 and 5.8 0.3, the best fit, based Eq. 15, Eqs. and 
experimental data achieved. Furthermore, does not vary function 
z/L. 

Table indicates that Eq. better matches experimental data than Eq. 12. 
Eq. subject some numerical error, but the linear density profile 


TABLE 3.—Statistical Evaluation Eqs. 12, and 


Control 
Compar- 


ison with 
Unstra- 
tified 
Comparisons Data 
with All Stratified and Equa- 


Data Having z/L 0.025 tion 


Statistical parameter 


(1) 


Standard deviation 


Calculated Student’s 
Tabulated Student’s 
Probability 

Number measurements, 


approximation buffers measurement errors affecting which turn affects 
Eq. 12. Since was based linear approximation, better curve-fitting 
techniques for the evaluation would improve the agreement between 
Eq. and experimental data. 

Finally, 5.8 0.3 comparable the accepted results Webb (20) 
for profiles the atmospheric boundary layer which 5.2 0.5. Previous 
evaluations range from (10). This variation results from two causes: 
(1) Taylor (17) notes, 0.6 and other small evaluations result from the 
consideration large and (2) and are not sensitive small changes 
the other. This insensitivity leads some divergence results. can also 
noted that Webb’s results follow from stratified flow studies which the 
scaling length, and buoyancy flux constant over the depth, whereas the 


studies presented here are characterized variable buoyancy flux over the 
depth. 


985 
(15) 

n 

(5) 

9.60 0.02 0.58 0.37 
1.65 1.65 1.65 1.67 
0.1 0.1 0.1 0.1 
109 109 109 
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The importance the width-to-depth ratio stratified flow noted from 
previous studies. Stratified flows narrow flumes are characterized peculiar 
vertical velocity profiles the width-to-depth ratio 4.3 less. addition, 
3-dimensional effects occur wider flows when the stratification highly stable. 

The studies presented here have demonstrated the application the Monin- 
Obukhov analysis predicting velocity gradients and velocity profiles stably 
stratified flows with 0.025. for paired data indicates the Prandtl-von 
Karman velocity law not applicable stratified flows, whereas velocity profile 
expressions, modified with the Monin-Obukhov scaling length, describe 
velocity gradients and vertical profiles stably stratified flows. Improved 
expressions for can lead better description vertical mixing. 

The coefficient for linear approximation the Monin-Obukhov function 
has been evaluated 5.8 0.3, which comparable the widely accepted 
value 5.2 0.5 for atmospheric boundary layer flows. 

Finally, these results demonstrate the importance stratified, open- 
channel flows and indicate that turbulent, 2-dimensional, stratified flows over 
rough boundary open channel are similar atmospheric boundary layer 
flows. Previous work French (4) narrow flume not valid for 2-dimensional 
stratified flows. French incorrectly defined the Monin-Obukhov scaling length. 
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buoyancy flux; 


ratio mass momentum mixing lengths; 

f(z) describing density gradient; 


von Karman’s constant 0.4); 
Monin-Obukhov scaling length; 
momentum mixing length; 


mass mixing length; 
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number measurements; 
gross Richardson number; 
temperature water; 

ambient temperature; 


u(z) velocity horizontal direction; 
shear velocity; 
average velocity flume center line; 
fluctuation vertical velocity; 
distance downstream; 
distance above bottom; 
Monin-Obukhov power series coefficient; 
constant density gradient; 
temperature difference over depth; 
roughness scale; 
dimensionless distance, z/L; 
kinematic viscosity; 
density, overbar denotes depth average, prime denotes density 
fluctuation; and 
function z/L. 
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TURBULENCE FOR FLOW 
OVER DREDGED CHANNELS 


INTRODUCTION 


The dredging channels alluvial material common feature many 
marine and river engineering schemes. Particular examples are the lengthening 
and deepening port approach channels, due the increasing size modern 
bulk carriers, and the construction dredged trenches across estuaries and 
rivers order bury pipelines. Unfortunately, many such dredged channels 
are alined angle the flow sediment, and this alinement can produce 
high siltation rates. Often, channels dredged cohesive material have steep 
sides, and high siltation rates are encouraged flow separation. 

Clearly, calculation infill rates dredged channels important the 
economic viability particular engineering schemes, and many calculation 
methods have been suggested, varying from simplistic approaches with bulk 
flow parameters (4), more recent ones, based detailed mixing processes 
(2,6,10). There are, however, difficulties with some the modern approaches, 
particularly for steep-sided channels. Approaches based mixing length tech- 
niques produce unrealistic mixing parameters for sediments (2), while more 
advanced methods assume spatially constant distribution dissipation parame- 
ters for internal energy (10). 

The present paper examines the use computer methods reproduce flow 
and mixing patterns steep-sided channel, dredged right angles the 
main direction the flow and sediment transport. 


Basic 


The governing turbulent mean equations for turbulent flow over dredged 
cut can written the form 
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Note.—Discussion open until January 1982. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication November 13, 
1979. This paper part the Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 107, No. HY8, August, 1981. ISSN 
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using the notation Fig. which flow velocities the and 
directions, respectively; the fluid density; the acceleration due 
gravity; the pressure the fluid; and equivalent turbulent 
shear stress the horizontal plane. Other stress gradients have been neglected 
comparison with the term. 

Integration Eq. with respect shows that the pressure within the 
fluid depends upon the density the fluid, and the mean and turbulent 
characteristics the flow field. However, most river and marine situations, 


FIG. Notation Model 


channel side slopes are such that only the density factor significant even 
when channel side slopes are steep. The pressure distribution with depth 
the present work thus taken hydrostatic first approximation. The 
computer methods presented can, however, used iteratively, required, 
include both density and flow field nonuniformity. 

Eqs. after the drop the vertical momentum term may rewritten 
the form 
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which effective eddy viscosity coefficient which must specified, 

Eqs. are solved. Previous work (2,10) has used both mixing-length 

and one-equation turbulence modeling specify v,, but these methods pose 

difficulties for separated flow situations. Both closure methods will, however, 

used the present work are illustrate their strengths and weaknesses. 
Mixing-Length Closure.—The eddy viscosity given 


which specified mixing length which varies with 
Turbulent Kinetic Energy Model.—The eddy viscosity given (9) 


which specified dissipation length which varies with and 
the turbulent kinetic energy per unit mass, given 


L(k) 


and 
(=) 


the preceding, numerical constant with value approx 0.3 (9). 

require specification the /,, parameter and suitable boundary 
conditions. Previous work (2,9,10) suggests that the /,,, distributions can 
expressed the form 


over the flow depth. Previous work (2,9) has also indicated various functional 
forms for but universal distribution exists. The present work thus 
uses typical distribution given 
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The integration Eq. from (water surface elevation) gives 
The turbulent shear stress approximated 
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for 
Boundary conditions are required at: 


the absence surface wind, implying zero shear stress 0). The zero 
gradient used (7) for the equation. 

The channel bed 0). nonslip condition used, for simplicity. 
This applies height equal the equivalent bed roughness 
height according logarithmic equation, Eq. 18. The computer method 
thus implicitly contains value related discretization. condition 
also used preference somewhat arbitrarily specified value (7). 

Upstream boundary. specified velocity profile just upstream the dredged 
channel section used. 

Downstream boundary. nonreflection condition used given 


Eqs. 4-11 are solved mixed procedure 
(3). splitting technique used time reduce the dimensionality the 
problem while one-dimensional Galerkin finite-element approach used over 
the flow depth minimize numerical errors and enable complex bed geometry 

The splitting technique (11) consists dividing the differential operators 
and into two parts 


inwhich 


Eqs. and can thus integrated numerically from time level nAt 
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leap-frog scheme used for successive computations and and allows 
the use values during the computation Model time steps 
are controlled the explicit steps, Eqs. 14a and 14c, and require that 


which the minimum space step used the computations. 
The water depth each flow cross section divided into constant number 


finite elements, shown Fig. The velocity and turbulent kinetic energy 
functions are approximated the series 


which the summation taken over the elements each cross section; 
matrix shape functions for element (i); and vector 


FIG. 2.—Discretization Flow Field 


nodal and values. Substitution Eq. into Eq. leads the appearance 
residuals, which are minimized the Galerkin approach, 


use linear-shape functions Eq. leads set equations for the 
unknown and values over the flow depth each cross section. 

Derivative terms appearing Eq. are replaced central finite-difference 
expressions, and the total set equations solved iteratively using Gauss-Seidel 
approach. Rapid convergence found, and computation times are kept low. 


steady uniform flow, possible derive analytical expressions for 
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and assuming velocity profile over the flow depth. logarithmic 
form used, e.g., according 


which where the bed shear stress; and von Karman’s 
constant; and and are related equations 


Analytical expressions for and are also known, once specified, 
namely 


the and, consequently, the analytical distributions, Eqs. 20-22, are 
based the drop the vertical diffusion term Eq. 9b, the distribution 


=H/20 size) 


computed 


u(cm/s) 


FIG. Roughness Height Implied Numerical Integration; 
Uniform Flow 


given Eq. will used this work. order use this distribution 
the computer solution, necessary minimize the difference between 
Eqs. 10, and and thus determine values and A,. One way 
achieving this minimization equating depth-mean differences zero, i.e., 
for the mixing length model 


u 
eight 
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K(z + H) 


which leads value 0.125. similar equation for the TKE model 
gives value 0.085. check the accuracy the proposed minimization 


860 


FIG. 4.—Real Flow Field and Its Discretization 


FIG. 5.—Determination Flume Roughness Height from Velocity Profile Uniform 
Flow 


method was provided operating the computer model with the calculated 
Values and comparing the velocity distribution obtained with Eq. 18. 
Fig. shows example this comparison. Use Eq. also enables 
determination the roughness height implied the computer model solved 
from the bottom the surface. Numerical tests with equally spaced 
elements showed that, for the mixing length closure, the implicit was found 
one-tenth the length the elements. value one-fifth the elements’ 
length was found for the closure. 
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FIG. 6.—Computed Longitudinal Variations 


comptd 
measrd 


FIG. Versus Computed Velocity Profiles; Mixing Length Closure 


FIG. 8.—Measured Versus Computed Velocity Profiles; TKE Closure 
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nonuniform flow, the A,, and factors vary spatially, and ideally, 
two-equation turbulence closure should used. However, first approxima- 


FIG. Versus Computed Eddy Viscosity Distributions; Mixing Length 
Closure 


FIG. 10.—Measured Versus Computed Eddy Viscosity Distributions; TKE Closure 
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FIG. 11.—Measured Versus Computed Shear Stresses 
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tion, the aforementioned approach can used, provided that adjusted 
locally. Integration Eq. over the flow mean depth leads the result 


Use Eqs. and then leads the result 


Ku, 


ait 


Similarly equation for can found 


Ku, 


The and values for use Eqs. 24-26 are obtained from the computer 
model and process successive approximations established determine 


FIG. 12.—Computed TKE Distributions 


the and values. Subsequent model tests showed that the method converged 
rapidly. 


order test the numerical procedures outlined previously, use has been 
made set laboratory experiments steady flow over steep trenches 
conducted the Delft Hydraulics Laboratory (8). The geometrical description 
the channel given Fig. Instantaneous velocities and shear stresses 
were measured various stations along the channel using laser doppler 
anemometer. The measured velocity profile station Fig. was used 
the upstream model boundary while the nonreflection condition Eq. 
was used the downstream boundary, station 11. The equivalent roughness 
height, z,, the channel was determined based the semilogarithmic 
plot the station velocity data (see Fig. 5). 


HY8 TURBULENCE MODEL 999 


Eqs. 24-26 were used with the computer velocities and produced 
variation through the channel, shown Fig. The velocity profile associated 
with Fig. for the mixing length closure shown Fig. from which 
can seen that the hydrostatic approximation represents good approach 
the steep-sided channel problem. The corresponding velocity figure for the 
closure shown Fig. the degree-of-fit both closure models seen 
very similar, and there little choose between them. 

Figs. and show the comparison between the two closure models for 
The disadvantage the /,, closure clearly seen stations and near 
the separated region, where very small values are predicted. The closure, 
the other hand, shows more realistic values the separated area and slightly 
higher values the accelerating part the flow; the large surface values are 
consequence the surface boundary condition. 

Reproduction the internal details the flow both methods rather 
similar, shown Fig. 11, for the shear stress term, The greatest area 
discrepancy occurs the downstream channel slope, due the development 
large values, shown Fig. 12. Generally, the closure model predicts 
larger t,, values than the /,, model, and has marginally better degree-of-fit 
the laboratory observations. 


The following conclusions may made: 


Steady flow across steep-sided dredged channels can predicted closely 
computer methods based fundamental flow equations, using hydrostatic 
pressure approximation. The hydrostatic pressure approximation works better 
stations and where the main flow nearly horizontal, than stations 
and where the flow converges, and there more pronounced vertical 
component the main flow. 

modeling both uniform and nonuniform flows reasonable cost and with 
reasonable accuracy. 

The computer method can readily used either with mixing length 
one-equation turbulence models, although the latter more advanced closure 
method requires considerably larger computation times. 

One-equation turbulence, closure provides better fit laboratory data 
than the mixing length approach, and obviates the problem singular values. 

The computer method with closure gives good reproduction 
and values, and could used for problems involving suspended loads, but 
provides less accurate representation internal parameters such 
The suggested method requires modification the postulation the 
use two-equation turbulence closure for problems involving bed load. 


Properties Base Functions 


(27) 
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Analytical Form Residuals Minimization Conditions 
element length, /,, the following condition must hold 


At At 


n+1/2 


dz; 


orinconcise form [A]{u}* {a} 


which 


Or, inconcise form 
inwhich 
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(33) 
(the repeated indices, imply summation over 
4, 
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similar approach the aforementioned leads the turbulent kinetic energy 
equations 


Typical elements the preceding matrices are 


The derivatives with respect and functions are replaced 
central finite differences staggered grid (the and values are computed 


the cross sections, and the values are computed midpoints between 
the two sections). The global matrices [A], [B], [C], {a}, {c}, 
{d} are formed superposition the matrix elements with the same indices 
(5) over cross section. The integrations Eqs. 33, 34, 37, 38, 41-44 are 
performed either analytically the aid Eq. numerically the use 
Gauss quadrature. 


Integration the Continuity Equation 
After the computation the profiles and their depth mean values, 
each cross section, the finite difference form Eq. becomes 
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The following symbols are used this paper: 


a 


turbulence closure parameter; 
normalized mixing length dissipation length distribution; 
total water depth; 

mean water depth; 

turbulence kinetic energy; 

dissipation length parameter; 

mixing length parameter; 

shape function; 

pressure; 

residual function; 

depth-mean velocity value; 

water velocity components; 

friction velocity; 

numerical integration steps (time, space); 
nondimensional water depth; 

Von Karman constant; 

mixing length, dissipation length factors; 
kinematic eddy viscosity; 

free surface elevation; 

shear stress turbulent flow; and 
bottom shear stress. 
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ROUTING DIVERGING- 
CONVERGING WATERSHEDS 


Necati 


INTRODUCTION 


Natural watersheds have complex geometric configurations. For modeling 
runoff natural watershed using physically based approach, necessary 
transform its complex geometry into simpler one having hydrologic response. 
The simplified geometric representation usually composed some combination 
three basic geometric elements: (1) Plane section (3,10); (2) converging section 
(13,15); and (3) diverging section (1,14). 

Some watersheds, both rural and urban, are roughly leaf-shaped. Their upland 
portions diverge, and their lower portions converge. watershed diverging 
and converging, then coupling the diverging section with the converging section 
may desired represent it. This geometric representation can called 
diverging-converging configuration. The questions are: does this representation 
allow better hydrograph prediction than alternative representations? And what 
are the effects varying geometric parameters the runoff hydrograph? 

this study, the overland flow equations depending the kinematic wave 
theory are solved using implicit finite difference method for the diverging- 
converging configuration. This diverging-converging flow model applied 
impervious surfaces which preserve the area and the length the watershed 
and fix the kinematic parameters for each event. The effects varying geometric 
parameters runoff hydrographs are investigated, and the results are compared 
with other geometric representations. 

The geometry diverging-converging surface shown Fig. which 
the length diverging section; the length flow for diverging 
section; the radial angle for diverging section; the degree divergence; 
the length converging section; the radial angle for converging 
section; the length flow for converging section; the degree 
convergence; the total flow length; and the total length the 
sections. The degree divergence, dimensionless parameter (R, 
L,)/R,. Also, the degree convergence, dimensionless parameter 
L,)/R,. 
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The basic equations unsteady flow open channel are the continuity 
and momentum equations. For diverging section, the continuity equation can 
written (14) 


which water depth; discharge per unit width; rainfall excess; 
time coordinate; and space coordinate. 


—v— 


FIG. 1.—Geometry Diverging-Converging Surface 


For converging surface, the continuity equation can expressed (13) 


However, for rectangular section, the continuity equation written 


The kinematic wave theory provides good approximation watershed flow 
due its large slope (7,10,11). For the kinematic approximation the momentum 
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equation can given 


which the friction slope; and the bed slope. Friction slope 
related flow regime using the Darcy-Weisbach, Chezy, Manning formulas. 
this study, the Darcy-Weisbach friction formula used, 


which the friction factor; the hydraulic radius; and the gravitational 
acceleration. 

Based upon some investigations (2,12), the friction factor can expressed 
for all practical flow regimes 


which constants depending surface characteristics; and the 
Reynolds number. Assuming for wide section, and putting g/v, 
elimination and from Eqs. 4-6 gives 


1/2-6 


which the kinematic viscosity. For conciseness, Eq. can written 


which and coefficients whose value depends the roughness 
the channel and the friction slope. Thus, for the kinematic wave approximation, 
the momentum equation becomes discharge-depth relationship. 

The substitution Eq. into Eq. gives 


This the kinematic wave equation for diverging flow. Similarly, Eq. yields 
the kinematic wave equation for converging flow 

calculate the flow from diverging-converging watershed, Eqs. and 
can orderly solved diverging and converging sections for boundary 
and initial conditions. For diverging section, the boundary condition 
when dR,. The assumption the downstream width diverging 
section equal the upstream width converging section results the boundary 
condition converging section equal the downstream condition diverging 
section, which, turn, results from the computation the first section. 
Therefore, for converging section, the upstream depth may equal the 
downstream depth diverging section. Initial conditions are the same for 
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For different geometric configurations, more useful compare the 
discharge values instead The discharge can expressed B,, 
which the discharge; and the width the watershed outlet. 
Therefore, the computation the discharge values requires the determination 
the outlet width which can estimated from other geometric parameters. 

The area the diverging section can written 


360 


which the area diverging section. The downstream width this 
section 


which the downstream width diverging section. Setting 


24, 


Assuming that the downstream width diverging section equal the 
upstream width converging section, the upstream width converging 
section can written 


watershed area; and Similarly, the outlet width can expressed 


2cA, 2cA 


2 


Because radial symmetry, and not enter explicitly into the 
computation. Thus, for preserving watershed area and length, there are only 
three geometric parameters, and which enter into the computation. 
These dimensionless parameters assume values between and depending 
watershed geometry. For some special values these parameters, the diverging- 
converging configuration gives some other idealistic geometric representations 
shown Table These six geometric configurations will represent most 
natural watersheds. 


Numericat 


Explicit and implicit finite difference methods have been used successfully 
many investigations kinematic wave modeling watershed runoff (1,4,6). 
Although explicit schemes are simple formulate and are easy program, 
they are conditionally stable. However, the implicit schemes are considerably 
more difficult program than explicit schemes, and they are theoretically 
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unconditionally stable linear sense. But they also have some practical 
limitations stability. 

Agiralioglu and Singh (1) used first-order and second-order explicit schemes 
for diverging and rectangular sections. But these schemes were not stable for 
converging sections, especially for small values degree convergence. Herein, 
four-point implicit scheme used for numerical solutions. Any dependent 


TABLE 1.—Geometric Parameters for Diverging-Converging Model and for Its Special 
Cases 


Geometric Parameters 


Degree Degree 
Geometric io, divergence, convergence, 
(1) 
Diverging-converging 
Diverging-rectangular 
Rectangular-converging 
Diverging 
Converging 
Rectangular 


O Known 
X Unknown 


FIG. 2.—Definition Sketch for Implicit Scheme 


variable, f(r,t), and its partial derivatives can approximated with four-point 
implicit finite differences follows (see Fig. 


which the space-weighting factor; and the time-weighting factor, 
respectively. 


j 
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Since the space coordinate, varies only space, the following expression 
holds for all values 


Similarly, assuming the rainfall excess only time-variant but spatially uniform, 
the following expression can written for all values 


Since the linear scheme, which good approximation overland flow 
(4), provides some simplification programming and saves computer time, 
linear routing method used. must noted, however, that the linear 
scheme requires more care selecting than nonlinear scheme. the 


linear routing method, the term can expressed the following finite difference 
form 


Substitution Eqs. 19, and into Eq. yields 


which the linear scheme equation for diverging section. 
Similarly, Eq. gives the following scheme for converging section 


+h*.,\° 


The kinematic diverging-converging model applied three events 
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impervious surfaces which are selected from the literature (5,7,9). The general 
characteristics each overland surface and its source are given Table 
which the width watershed; and the rainfall duration. 

The kinematic parameters, and are estimated from the reported data. 
For example, Morgali specified 1.0 and 36.3 for the rising hydrograph 


TABLE 2.—Characteristics Selected Events 


Reference 
(2) 
Morgali (1970) 
Langford and 

Turner 
(1973) 
Muzik (1973) 


Assumed. 


TABLE 3.—Estimated Kinematic Parameters 


Kine- 

matic 

vis- 
cosity, 

millions 

Parameters for Parameters for 
meters 

per Rising Hydrograph Recession Hydrograph 


second 
(2) (7) 
1.30 36.30 1.00 16,630.64 32.00 1.00 
1.14 5.40 0.24 2.16 1.70 76.60 0.75 
1.00 32.00 1.00 509,899.28 3.00 24.00 1.00 


and 32.0 for the recession hydrograph. this case, and 16630.64 
for the rising hydrograph and 18865.38 for the recession hydrograph for 
1.3 m’/s. These kinematic parameters are indicated Table 
All computations depend the assumption constant kinematic parameters 
for different geometric representations. This assumption may reasonable 


Rain- 
fall Rain- 
den- fall 
sity, 
ginal Width Length meters 
Serial test Paved B,in Slope, per sec- 
numbe number type meters| meters hour onds 
(1) (3) (4) (5) (6) (7) (8) (9) 
sheet 
Se- 
rial 
num- 
ber 
(1) 
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because the kinematic parameters are related flow regime but not geometric 
parameters, except bottom slope seen from and The bottom slope 
considered constant for each surface. 

The weighting values, 0.5 and 0.5, were used for all computations. 
For the accuracy the results, should small, and for the achievement 


: RECTANGULAR 
(MEASURED) 


t RECTANGULAR 
(COMPUTED) 


: DIVERGING-CON- 
VERGING (O=05, 
c=0.1, d=00) 


DISCHARGE, 


3 4 
TIME, t, MINUTES 


FIG. 3.—Comparison Measured and Computed Hydrographs 


: DIVERGING- CONVERGING 
( p=05, d=0.0, c=01) 


: DIVERGING 
( p=! 0, d=01 ) 


CONVERGING 


OISCHARGE, 


TIME, t, MINUTES 


FIG. 4.—Comparison Hydrographs for Different Geometric Representations 


the stability practice, should smaller than certain limit. Therefore, 
the space interval was selected one-fortieth the flow length, L/40. 
The time interval was selected, and set 2.0 sec, 0.5 sec, and 0.2 sec for 
Langford and Turner’s, Morgali’s, and Muzik’s data, respectively. Larger values 
have been tested, but they have given results which deviated from the 
expected values for equilibrium condition. 


4 
$00 
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200 \ 
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2000; 
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The hydrograph obtained from Morgali’s data for diverging-converging 
configuration for 0.5, 0.0, and 0.1 compared with the observed 
hydrograph and the computed hydrograph for rectangular section shown 


DISCHARGE, 


40 so 
TIME, ¢, SECONDS 


FIG. 5.—Effect Degree Divergence Equilibrium Hydrograph for Diverging- 
Converging Model 


300 


° 
(2) 


DISCHARGE, 


TIME, MINUTES 


FIG. 6.—Effect Degree Convergence Partial Equilibrium Hydrograph for 
Diverging-Converging Model 


Fig. There good agreement between the equilibrium measured and 
computed discharges. This indicates that the implicit scheme with linear routing 
good approximation for the computations. 
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Fig. shows the effect watershed geometry the runoff hydrograph 
with Langford and Turner’s figures for diverging-converging, diverging, and 
converging configurations. seen from this figure, the rising hydrograph for 
diverging configuration goes more rapidly the beginning than the rising 
hydrograph for converging section. the contrary, the recession hydrograph 
goes down more rapidly for diverging configuration than the recession 
hydrograph for converging configuration. The runoff hydrograph for diverg- 


OIMENSIONLESS DISCHARGE 


OIMENSIONLESS TIME 


FIG. 7.—Effect Degree Divergence Dimensionless Rising Hydrograph 


DIMENSIONLESS D0! SCHARGE 


04 06 08 
OIMENSIONLESS TIME 


FIG. 8.—Effect Degree Convergence Dimensionless Rising Hydrograph 


ing-converging configuration stays between the values these two configurations. 
the other hand, the diverging representation gives the longest time 
concentration. 

Fig. shows the effect degree divergence, the runoff hydrograph 
for Muzik’s data. this case, the hydrographs not change significantly 
the value increases. 
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obtain partial equilibrium hydrograph, the rainfall duration, was selected 
sec for Morgali’s data. For this case, the effect the degree convergence, 
the runoff hydrograph given Fig. This figure shows that the 
peak discharge varies with geometric parameters the partial equilibrium 
condition. 

general, the results are found good agreement with previous 
study, which diverging and rectangular flows were investigated impervious 
surfaces using explicit finite different methods (1). 


: CONVERGING (c=0.1) 


: DIVERGING-~ CONVER- 
GING 


DISCHARGE 


TIME 


FIG. 9.—Comparison Dimensionless Recession Hydrographs for Different Water- 
shed Geometries 


order show clearly the effect watershed geometry both the rising 
and recession curves, dimensionless discharge and time variables are designated 
with asterisks 


25) 


which the normalizing discharge which equals the equilibrium discharge, 
iA; and the normalizing time which equals the time concentration 
for rectangular section, i.e., 

Fig. shows the dimensionless rising curves for different values including 
0.0, 0.5, and 0.9 for 0.5 and 1.0. This figure clearly 
the effect the rising hydrograph than Fig. 
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Fig. shows number dimensionless rising hydrographs from Morgali’s 
for different values 0.1, 0.5, and 0.9, for 0.0 and 0.5. 
diverging-converging configuration, the parameter efficient influencing 
the shape hydrograph, but has little influence the time concentration. 

Fig. shows the effect the parameter, the recession hydrograph 
from Muzik’s data. seen from this figure, the recession hydrographs 
down more rapidly for larger values than recession hydrographs for smaller 
values The diverging configuration 1.0) can give smaller lag time 
than the others. Fig. supports this conclusion, also. 

Fig. indicates the rising hydrograph from Morgali’s data for different 
geometric representations, viz., diverging-converging, converging, diverging, and 
rectangular representations. examination this figure shows that the 
converging representation gives the shortest time concentration, while the 


: CONVERGING (c=01) 
: DIVERGING (d=00) 
: RECTANGULAR 

DIVERGING-CONVER- 


GING P=05, 
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FIG. 10.—Comparison Dimensionless Rising Hydrographs for Different Watershed 
Geometries 


diverging representation gives the longest one. The time concentration for 


the rectangular representation stays between the values these two representa- 
tions. 


Based the numerical study kinematic flow equations, the following 
conclusions are drawn: 


diverging-converging surface, the geometric parameters are important 
factors both rising and recession hydrographs. diverging-converging flow 
model gives hydrograph values between values diverging and converging 
model depending geometric parameters. For some special values geometric 
parameters, diverging-converging models can used for other idealistic geometric 
representations. 

Different geometric representations give runoff hydrographs with different 
shapes under the same conditions. comparison three basic geometric 
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representations (diverging, rectangular, and converging) shows that diverging 
representation gives the longest time concentration, while converging 
representation yields the shortest one. The reverse true for lag time. 

simulate the hydrographs for leaf-shaped watershed, diverging-con- 
verging representation may better than the other idealistic representations. 

Since the stage, changes more slowly for diverging representation 
than others, diverging representation the most stable. the other hand, 
converging representation the least stable. Computation converging 
section needs small time intervals for stability, especially for small values 
the degree convergence. However, the four-point implicit method can 
used successfully for various geometric representations watershed with 
large variations stage. 
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The following symbols are used this paper: 


watershed area diverging section; 
watershed area converging section; 
constant; 
average width watershed; 
downstream width diverging section; 
downstream width converging section; 
constant; 
degree convergence for converging section; 
degree divergence for diverging section; 
Darcy-Weisbach friction factor; 
acceleration due gravity; 
water depth; 
rainfall excess; 
space discretization index; 
time discretization index; 
length flow; 
length diverging part region; 
length converging part region; 
discharge; 
normalizing discharge; 
dimensionless discharge; 
discharge per unit width; 
Reynolds number; 
length diverging section; 
length converging section; 
hydraulic radius; 
space coordinate; 
slope friction; 
slope bed; 
rainfall duration; 
normalizing time; 
dimensionless time; 
time; 
weighting factor spatial derivative; 
weighting factor temporal derivative; 
total length sections; 
routing coefficient; 
routing coefficient; 
spatial increment; 
temporal increment; 


d= 
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radial angle for diverging section; 
radial angle for converging section; 
kinematic viscosity water; and 
length ratio. 
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DISCUSSION 


Note.—This paper part the Journal the Hydraulics Division, Proceedings 
the American Society Civil Engineers, Vol. 107, No. HY8, August, 1981. 
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Frequency 


Discussion Arun Kumar’ and Rema 


The author has proposed Bayesian regression analysis for predicting flood 
order account for uncertainties model selection and the evaluation 
parameters from limited data. The method illustrated fitting two-parameter 
distributions (normal, log normal, and Gumbel) observed flood data. is, 
however, well-recognized that three-parameter distributions such Pearson 
Type III, log Pearson Type III, and general extreme distribution functions fit 
better flood data than two-parameter distributions (5,13). the writers’ 
opinion, the proposed methodology has very limited applicability when the 
three-parameter distribution fitted. Unlike two-parameter distribution, con- 
struction probability paper for three-parameter distribution will depend 
upon the coefficient skewness observed sample. other words, for 
given distribution, sampling errors the estimation skewness will result 
different probability scale papers. Even additional data may cause change 
the probability scale. Since the proposed method does not account for 
change probability scale for given distribution, uncertainty resulting 
the estimation skewness due sampling errors cannot accounted for. 

Gumbel, normal, and log normal distributions have been fitted Salt Fork 
data and has been pointed out that although the log normal distribution apparently 
shows the best fit probability paper, may not appropriate distribution 
for predicting flood large return periods. The writers feel that the conclusion 
derived the author not correct. 

The goodness fit different distributions can checked comparing 
the observed and predicted values various return periods. The data considered 
follow that distribution which gives the closest fit (5,12). The comparative 
goodness fit Salt Fork data using the normal, the log normal, and the 
Gumbel distribution given Table observed from the table that 
the log normal distribution gives better fit compared the normal and 
the Gumbel distribution, because has the least-mean square error and largest 
correlation coefficient. 

The log normal distribution assumes apriori that the logarithm flows follow 
normal distribution. This may not true all cases. such cases, the 
goodness fit can further improved using data-based power transforma- 
tion (11), which reconstitutes the data such that the transformed series follows 
normal distribution. This transformation which given 


1980, Wilson Tang (Proc. Paper 15532). 


Prof. Civ. Engrg., Delhi Coll. Engrg., Kashmere Gate, Delhi, 110006, 
India. 


Asst. Prof. Civ. Engrg., Indian Inst. Tech., New Delhi, India. 
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transforms the original series, Y,, into approximately normal series, Z,, 
for some value can seen from Eq. that log normal special 
case power transformation. For Salt Fork data, when —0.14, the mean, 
standard deviation, skewness, and kurtosis the transformed series, Z,, are 
4.7533941, 0.1769088, —0.0015, and 3.5, respectively. Taking the transformed 
series normal, the flood peaks various return periods are calculated and 
given Table observed that the power transformation gives even 
better fit than does the log normal distribution. The 20-yr flood using the power 
transformation 7,047 with standard deviation 1,155. 


TABLE 6.—Observed and Fitted Flood Peaks Salt Fork (17 yr) Using Various Models, 
cubic feet per second 


Fitted Flows 


Observed flows Log normal Power transformation 
(1) (4) (5) 


Note: For 1/(n the values for the various fitted 
flows are 683, 496, 428, and 254, respectively. For the correlation coefficient, the values 
are 0.911, 0.954, 0.973, 0.977, respectively. 


The author seems infer that the normal distribution, which has less 
standard deviation should preferred comparison the log normal and 
the power transformation for the estimation flood large return period. 
The standard deviation, given Eq. 18, depends upon the type distribution 
chosen. can index for selecting distribution only when all the competing 
distributions fit the data equally well. the writers’ opinion, since the variation 
the Salt Fork data more accurately expressed the log normal and the 


power transformation, the extrapolation flood should estimated using 
them. 
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6,860 5,660 6,382 6,507 6,806 
6,030 5,016 5,222 5,312 
5,380 4,579 4,665 4,497 4,510 
3,370 4,229 4,185 3,989 3,967 
2,990 3,927 3,807 3,597 3,556 
2,750 3,653 3,481 3,275 3,226 
2,640 3,397 3,193 3,000 2,949 
2,630 3,151 2,931 2,758 2,708 
2,600 2,911 2,687 2,541 2,494 
2,350 2,670 2,454 2,340 2,299 
2,340 2,425 2,229 2,151 2,117 
2,300 2,168 2,005 1,971 1,945 
1,920 1,895 1,777 1,794 1,778 
1,800 1,592 1,540 1,618 1,613 
1,350 1,242 1,280 1,435 1,443 
1,230 805 978 1,236 1,258 
950 162 572 992 1,034 
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Chander, S., Spolia, K., and Kumar, A., Frequency Analysis Power 
Journal the Hydraulics Division, ASCE, Vol. 104, No. HY11, 
Proc. Paper 14179, Nov., 1978, pp. 1495-1504. 


12. Gupta, L., Frequency Distribution Water Resources Research, 
Vol. 1970, pp. 1193-1198. 


13. Jenkinson, F., Maximum Technical Note 98, World 
Meteorological Office, Chapter 1969, pp. 183-257. 


The author presents new method fitting hydrologic data linear regression 
probability distributions which are plotted straight lines specific 
probability papers. This method permits combining flood estimates from different 
distributions under Bayesian framework. The writer offers the following 
comments. 

Because the proposed fitting procedure innovative, very little known 
about its statistical properties. Thus, before advocating the method effective 
fitting procedure, there need to: (1) Establish its statistical properties, such 
its bias and efficiency estimating different design events; and (2) compare 
these properties with the traditional fitting methods moments (MO) and 
maximum likelihood (ML). Table the writer presents the 100-yr design 
values for the author’s sample and ML. These results show that the 
author’s method overestimates the design values compared and 
methods, especially for the log normal and Gumbel models, where the differences 
are high 9%-25%. present, there way telling whether the 
author’s results are improvement over the and methods. However, 
since the properties and are relatively well-known, the high disparities 
the present results suggest, best, that probably the author’s method will 
produce biased (and thus, poor) estimates design values for most models. 

The author makes understatement saying ‘‘the Monte Carlo simulation 
studies assume the existence known probability distribution the real 
Monte Carlo simulation powerful technique available with the advent 
digital computers systematically study several properties the data samples 
and fitting methods random variable; such properties are ordinarily not 
tractable mathematically. This technique, e.g., the only means available 
study the properties the author’s proposed fitting method. With regard 
the relevance the Monte Carlo simulation real world situation, samples 
from several theoretic distributions exhibit similar properties, such that the 
sample skew small size sample underestimate the population skew 
(19) that the estimates design values are less variable than other 
methods, reasonable believe that the samples from the unknown distribution 
the real world also will behave accordingly. 

Finally, the concept that each model could potentially correct may not 
acceptable this time. Historically, many methods were evolved and used 


1980, Wilson Tang (Proc. Paper 15532). 
St. Johns River Water Management Dist., P.O. Box 1429, Palatka, Fla. 32077. 
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derive flood magnitude-frequency relations (15,17). the course time, 
several methods were discontinued and others had varying degrees popular 
acceptance. possible reason why each several distributions had strong 
proponents among hydrologists the past, that humid area where flows 
are less variable, all methods gave results that are generally satisfactory (16). 
Also, some agencies chose particular model and continued use either 
for consistency due certain theoretic appeal the models (i.e., Gumbel 
the United States Weather Bureau, now the National Oceanic and Atmospheric 
Administration and log normal the United States Army Corps Engineers). 
Eventually, more data became available from many stream gaging stations, 
studies were conducted compare the results given different models 
general fashion (14,16); such studies indicate that some models give inaccurate 
(i.e., unsatisfactory) results under certain conditions. For given hydrologic 
sample, reasoning, several models could eliminated unfit follows. 
Two basic attributes probability model are: (1) The distribution 
probabilities the random variable, i.e., how the cumulative distribution function 
plots reference probability paper, like normal log normal; and (2) the 
statistical parameters, SP, the random variable, such mean, variance, 


TABLE 7.—Predicted Design Values Methods Moments and Maximum Likeli- 
hood, cubic feet per second 


Data from 1959-1975 Data from 1959-1977 


method method method method 
(2) (3) (4) (5) 


Note: cfs 0.028 m’/s. 


Normal 
Log Normal 
Gumbel 


and the coefficient skewness, CS. For the real world data, the cumulative 
distribution function not known; however the estimates can obtained 
from sample data. Even though these estimates could biased, especially 
that CS, they still provide qualitative information about the real world random 
variable. For the 1959-1975 and 1959-1977 data Salt Fork (author’s sample), 
the conventional unbiased estimates skewness are: 1.354 and 1.511, 
respectively, for data real space, and 0.193 and 0.379, respectively, 
for data log space. Since the size the samples small, these values 
may regarded underestimates population values. the basis these 
estimates, none the models selected the author valid; for validity the 
normal distribution requires real space, the log normal distribution 
requires log space, and the Gumbel distribution requires 1.14 
real space. The design values predicted the three models differed greatly 
the and methods (see Table 7), well the author’s fitting 
method (see Tables and the original paper) because erroneous models 
were selected. The models are three-parameter probability distributions 
(18). 

described Ref. 18, the method (using the values real 


Model 

(1) 

6,780 6,780 6,540 6,540 

8,700 8,440 8,410 7,860 
8,130 7,250 7,840 6,820 
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space mentioned the foregoing and other SP), the writer evaluated the 100-yr 
design values four three-parameter models, i.e., log normal, Weibull, Pearson 
type III, and log Pearson type III. The results obtained these four models 
are: (1) 8,290 cfs, 8,290 cfs, 8,300 cfs, and 8,320 cfs, respectively for 1959-1975 
data cfs 0.028 and (2) 8,130 cfs, 8,160 cfs, 8,160 cfs, and 8,160 
cfs, respectively for 1959-1977 data. Thus, the four models produce practically 
equivalent results. Probably, the author’s fitting method will also lead similar 
results. However, the author’s method rather tedious since for each value 
separate probability paper has designed the case three-parameter 
model. The writer’s choice the four models shows that under appropriate 
choice models, the model uncertainty greatly reduced. 

summary, the author’s fitting procedure needs investigation for its 
statistical estimation properties. The severe model uncertainty portrayed the 
author the result erroneous choice models; under appropriate 
choice models all models might produce practically equivalent results, 
which case the Bayesian approach ineffective its intended purpose. 
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Errata 


The following corrections should made the original paper: 


Page 1366, Eq. 2b: Should read instead 


2m +t 2ms+t 
1980, Guus Stelling (Proc. Paper 15660). 
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GEOMETRY GRAVEL 


The author makes the reasonable assumption that self-adjustment gravel- 
bed stream takes place response independently-determined values water 
flow rate, and sediment transport rate, G,. The sediment transport rate 
then related the bed shear stress, the channel. 

Bagnold (24,25) has shown that correlates extremely well with stream 
power QS. The author has himself used this relationship previously (3). 

thus closely related G,, which independently determined, 
clear that not free minimized the author hypothesizes. 
The writer would appreciate the author’s comments this point, which 
relevant much the author’s previous work (3,4,5). 


24. Bagnold, A., Approach the Sediment Transport Problem from General 
Professional Paper United States Geological Survey, Washington, 
D.C., 1966. 


Bagnold, A., Transport Natural Rivers,’’ Water Resources Research, 
Vol. 13, No. Apr., 1977, pp. 303-312. 


Minimum Power: 


The quantitative credibility the minimum stream-power hypothesis depends 
part the second author’s (16) use measured data support his proposition 
that unit stream power, Vs, minimized alluvial channels. the present 
paper the authors claim that this special case minimizing yQs. Guy, 
Rathbun, and Richardson (23) and Maddock (25) infer that the flume data used 


“September, 1980, Howard Chang (Proc. Paper 15678). 


Lect., Dept. Agr. Engrg., Lincoln Coll., Univ. Canterbury, Canterbury, 
New Zealand. 


1980, Charles Song and Chih Ted Yang (Proc. Paper 15691). 
Lect., Dept. Agr. Engrg., Lincoln Coll., Univ. Canterbury, Canterbury, 
New Zealand. 


“Reader, Dept. Civ. Engrg., Univ. Canterbury, Christchurch, New Zealand. 
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Yang (16) support his proposition was acquired under conditions which 
either the discharge, and the water depth, were preselected and thus 
independent variables, and the energy slope, were preselected. Sediment 
and fluid properties were also constant within given set. The river data used 
presumably described short-term (periods days weeks) variations under 
conditions constant and 

Under each these two sets conditions, constraints, the hypothesis 
minimum seems questionable. and are independent variables, 
then the value imposed externally and cannot minimized adjustment 
the dependent variables. With and independent, Raichlen and Kennedy 
(28) and Mantz (26,27) have reported progress towards equilibrium laboratory 
channels with initially-plane beds. each case, the bed shear stress, and thus, 
energy slope and stream power, increased monotonically equilibrium was 
approached. Thus, achieving its equilibrium value stream power, the system 
passed through states lower stream power, each which the imposed 
constraints and were satisfied. does not seem reasonable that channel 
development ceased result stream power minimization. 

third set constraints may that and the sediment transport rate, 
are imposed the system, the long-term development natural 
streams. this case, minimization yQs implies minimization since 
imposed. single-valued function shear stress (29), then yds 
fixed, and minimization implies that maximum. The equilibrium 
would occur, under the minimum stream power hypothesis, maximum value 
(19,20), then with fixed, the stream power fixed and minimization 
possible. possible variation stream width, does not affect this argument, 
since function (24), and thus not variable when independent. 

alternative hypothesis proposed the writers (22) that self-formed 
stream achieves equilibrium when its resistance coefficient given 


local maximum. Examples from variety alluvial channel situations 
were cited support this view; Ref. also lends weight this. 

Considering the third set constraints with and independent, maximizing 
implies maximizing depends shear stress, then given and 
thus, maximized and minimized. This agreement with minimizing 
under these constraints. If, however, depends yQs, then fixed 
and maximizing implies maximizing thus maximizing gives reasonable 
result here, whereas minimization yQs could not occur. 

With and independent, maximizing implies maximizing and bed shear 
stress. This compatible with the increases shear stress reported (26,27,28) 
flume experiments approach equilibrium, whereas the application minimum 
stream power this case requires explanation. 

With and independent, maximizing implies maximizing This agrees 
with the behavior lower-regime bed forms (21), whereas yQs cannot 
minimized under these constraints. 


| 
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The hypothesis maximum thus appears more widely applicable, 
and some cases more realistic than that minimum 
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and descriptions must placed text rather than within the figure. 


Tables should typed (an original ribbon copy and two duplicates) one side 
(220-mm) (280-mm) paper. explanation each table must appear the text. 


References cited text should arranged alphabetical order appendix the 
end the paper, preceding the Appendix.—Notation, Appendix.—References. 


10. list key words and information retrieval abstract 175 words should provided 
with each paper. 


summary approximately words must accompany the paper. 
12. set conclusions must end the paper. 


13. Dual units, i.e., U.S. Customary followed (International System) units parentheses, 
should used throughout the paper. 


14. practical applications section should included also, appropriate. 
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